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ABSTRACT 

The intergalactic magnetic field (IGMF) may leave an imprint on the angular 
anisotropy of the extragalactic gamma-ray background through its effect on electro- 
magnetic cascades triggered by interactions between very high energy photons and the 
extragalactic background light. A strong IGMF will deflect secondary particles pro- 
duced in these cascades and will thus tend to isotropize lower energy cascade photons, 
thereby inducing a modulation in the anisotropy energy spectrum of the gamma-ray 
background. Here we present a simple, proof-of-concept calculation of the magnitude 
of this effect and demonstrate that current Fermi data already seem to prefer non- 
negligible IGMF values. The anisotropy energy spectrum of the Fermi gamma-ray 
background could thus be used as a probe of the IGMF strength. 

Key words: magnetic fields - diffuse radiation - gamma rays: theory - galaxies: 
active - BL Lacertae objects: general - galaxies: starburst 
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1 INTRODUCTION 

Cosmic magnetic fields are expected to play a fundamental 
role in the physics of a large variety of astrophysical systems. 
The fate of high energy particles, such as ultra-high energy 
cosmic rays (UHECRs), as they propagate through the Uni- 
verse and the acceleration of charged particles in astrophys- 
ical objects hinge, in large part, on magnetic fields within 
their sources and those permeating the Universe. Large-scale 
magnetic fields such as those found in galaxies and clus- 
ters of galaxies have been a mystery for the past several 
decades. While they are generally thought to be the result 
of the amplification of weak seed fields, the origins of these 
seed fields, whether they be cosmological or astrophysical 
in nature, are largely unknown. A definitive measurement 
of the intergalactic magnetic field (IGMF) could provide a 
fundamental step in resolving the questions of the origins of 
cosmic magnetic fields and their impact on the evolution of 
the systems in which they reside, but sufficiently constrain- 
ing observations have thus far remained elusive. Previously, 
upper bounds of ~ 10 -9 G have been found through Fara- 
day rotation limits of polarized radio emission from distant 
quasars and the study of the effect of magnetic fields on 
the anisotropy of the cosmic microwave background (for a 
brief review of observational and theoretical techniques to 
constrain the IGMF, see e.g., Neronov & Semikoz 2009, and 



references therein). The recent availability of data from the 
Fermi Large Area Telescope (LAT) has renewed interest in 
the IGMF as observations of the gamma-ray sky and its par- 
ticipating source populations provide a unique perspective 
on the strength and structure of the IGMF. 

The gamma-ray sky consists of resolved point sources 
(such as normal and active galaxies, pulsars, etc.), transient 
gamma-ray sources (e.g., gamma-ray bursts), and the diffuse 
gamma-ray radiation comprised of emission from the Galaxy 
and the isotropic (presumably, extragalactic) gamma-ray 
background (EGB). The EGB is a window into the high- 
energy processes in the Universe, and its origins have been 
the subject of much debate. It is expected that emission 
arising from unresolved, extragalactic point sources, such as 
blazars and starforming galaxies, comprises a sizable contri- 
bution to the EGB (see e.g., Stecker & Salamon 1996; Fields 
et al. 2010; Abdo et al. 2010d; Stecker & Venters 2011, and 
references therein). Additionally, many of these extragalac- 
tic point sources are also sources of very-high energy (VHE) 
gamma rays, 1 which interact with the soft photons of the 
extragalactic background light (EBL), consisting of the cos- 
mic microwave background (CMB) and the background of 
infrared, optical, and ultraviolet radiation from direct and 
dust-reprocessed starlight. The interactions of VHE pho- 
tons with the EBL initiate electromagnetic (EM) cascades, 
giving rise to another contribution to the EGB (see e.g., 
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In this paper, we take VHE to be ~ TeV. 
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Coppi & Aharonian 1997; Kneiske & Mannheim 2008; In- 
oue & Totani 2009; Venters 2010). However, the degree to 
which astrophysical sources contribute to the EGB depends 
on their distribution with respect to luminosity and redshift 
(the gamma-ray luminosity function, GLF), their intrinsic 
spectra at gamma-ray energies, and the nature of the EBL, 
all of which remain the subjects of intense debate (see e.g., 
Venters 2010; Fields et al. 2010; Abdo et al. 2010b,d; Makiya 
et al. 2011; Inoue 2011; Stecker & Venters 2011; Abazajian 
et al. 2011). 

A promising complementary technique for investigat- 
ing the contributions of astrophysical sources to the EGB 
involves studying the angular fluctuations in intensity of 
the EGB as a function of energy, the anisotropy en- 
ergy spectrum. As pointed out in Siegal-Gaskins & Pavli- 
dou (2009; see also Hensley et al. 2010; Siegal-Gaskins et al. 
2011), the contribution of a given class of sources to the 
overall fluctuation angular power at a given angular scale 
of the EGB at a given energy is weighted by the square 
of its fractional contribution to the intensity of the EGB 
at that energy. As such, if the relative contributions of the 
various EGB contributors change with energy, so too will 
the fluctuation angular power as a function of energy. Ex- 
pectations for the anisotropies of the commonly invoked 
astrophysical contributors to the EGB are that starform- 
ing galaxies contribute very little anisotropy owing to the 
fact that they are so numerous and individually faint, while 
blazars, being far more rare than starforming galaxies but 
much brighter, could have a large contribution to the EGB 
anisotropy (Ando et al. 2007; Ando & Pavlidou 2009; Siegal- 
Gaskins & Pavlidou 2009; Hensley et al. 2010; Ackermann 
et al. 2012a). It has also been suggested that annihilation 
in dark matter substructures in the Galactic halo could also 
have a large contribution to the EGB anisotropy, even if the 
contribution to the EGB intensity from dark matter anni- 
hilation is subdominant with respect to other contributions 
(see e.g., Ando et al. 2007; Siegal-Gaskins 2008; Taoso et al. 
2009; Ando 2009; Siegal-Gaskins & Pavlidou 2009; Hensley 
et al. 2010). Early anisotropy results from the Fermi-L AT 
Collaboration (Ackermann et al. 2012a) suggest that the 
energy dependence of the fluctuation angular power of the 
EGB may be consistent with arising from a source popula- 
tion with spectra similar to that of blazars (a ~ 2.4), even 
though blazars may comprise < 30% of the EGB intensity 
(Abdo et al. 2010d; Cuoco et al. 2012). Thus, even contribu- 
tions to the EGB that are subdominant with respect to other 
contributions can have a profound impact on the anisotropy 
of the EGB, especially if the relative prominence of its com- 
ponents changes with energy. As such, the anisotropy energy 
spectrum could prove a powerful tool in detecting changes in 
the EGB anisotropy resulting from the magnetic deflection 
of EM cascades from blazars. 

In interacting with EBL photons, VHE gamma rays 
initiate EM cascades of photons, electrons, and positrons. 
The interaction of a VHE gamma ray with an EBL photon 
produces a pair of electrons and positrons, which will, in 
turn, Inverse Compton scatter EBL photons to high ener- 
gies. These up-scattered photons will, in turn, pair produce, 
and the process continues until the energies of the resulting 
photons are low enough that pair production is no longer 
efficient. In the presence of an IGMF, the charged parti- 
cles in the cascade are deflected away from the direction of 



propagation of the primary photon. As these charged par- 
ticles propagate through the EBL, they can up-scatter soft 
photons towards the observer, though they will appear to be 
coming from a different direction from the GeV gamma rays. 
Observationally, this process results in a "halo" of gamma 
rays around the source. The detection of a halo around a 
gamma-ray point source would not only provide strong evi- 
dence for an IGMF, but would also provide an indication of 
the strength of the IGMF. Several groups have already con- 
ducted searches for extended emission around gamma-ray 
sources (see e.g., Aharonian et al. 2001; Neronov & Vovk 
2010; Ando & Kusenko 2010; Aleksic et al. 2010; Neronov 
et al. 2011; Arlen et al. 2012), though as yet, no significant 
indication of gamma-ray halos distinguishable from instru- 
mental effects have been reported. However, the absence of 
significant halo emission could imply a strong IGMF that 
would extend the EM cascade halo such that the surface 
brightness of the extended source falls below the detection 
threshold (Neronov & Vovk 2010; Ando & Kusenko 2010). 
In addition to possibly creating a halo, the deflection of EM 
cascades by the IGMF could result in a reduction in ob- 
servable cascade radiation. Analyses employing simultane- 
ous observations of blazars by atmospheric Cherenkov tele- 
scopes and Fermi make use of this fact to derive limits on the 
IGMF (Dermer et al. 2011; Essey et al. 2011; Taylor et al. 
2011; Kachelriess et al. 2011). These techniques are quite 
promising and, with more data from Fermi, will continue to 
shed light on the IGMF. In this paper, we consider a dif- 
ferent promising technique that makes use of the anisotropy 
energy spectrum of the EGB to constrain the IGMF. 

For any cosmological population, such as blazars, that 
emits gamma rays at very high energies, the effect of EM 
cascading results in a flux suppression at the highest en- 
ergies and enhancement at lower energies as the cascades 
redistribute radiation from the higher energies to the lower 
energies (see e.g., Strong et al. 1973, 1974; Coppi & Aha- 
ronian 1997; Kneiske & Mannheim 2008; Inoue & Totani 
2009; Venters 2010). In considering the spectral subpopula- 
tions of blazars, fiat-spectrum radio quasars (FSRQs) and 
BL-Lacertae type objects (BL Lacs), the intrinsic spectra 
of FSRQs are much softer than that of both BL Lacs and 
EM cascades, and their relative prominence with respect to 
the EGB can change with energy (Venters 2010; Venters 
& Pavlidou 2011). As such, the fluctuation angular power 
in the EGB should change with energy even if FSRQs and 
BL Lacs have similar anisotropy properties. Furthermore, as 
noted earlier, the deflection of EM cascades by the IGMF 
can spatially extend the emission around gamma-ray sources 
and reduce the observable cascade radiation. Thus, the im- 
pact of a strong IGMF is to reduce the blazar anisotropy 
both because of the increased angular scale of the sources 
and because of the reduction in the overall blazar contri- 
bution to the EGB (assuming that blazars are the most 
anisotropic EGB component). On the other hand, in the 
absence of a strong IGMF, the EM cascades are highly col- 
limated in the direction of the primary photon (angular 
spread ~ m e /E e ), and their impact on the blazar anisotropy 
is to reinforce it due to the enhancement in the overall blazar 
contribution to the EGB. Hence, in the two extreme cases 
for the IGMF, EM cascades could have a profound and oppo- 
site impact on the blazar anisotropy and, by extension, the 
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EGB anisotropy. Therefore, studies of the EGB anisotropy 
could prove a sensitive probe into the nature of the IGMF. 

In this paper, we explore the impact of EM cascade ra- 
diation from blazars on the EGB anisotropy in the extreme 
cases of a negligible IGMF and an IGMF that is sufficiently 
strong to essentially isotropize the cascade emission (i.e., 
from the observationally practical standpoint). This study 
is meant to serve as a proof of concept for the feasibility of 
using EGB anisotropics as an IGMF probe. In Section 2, we 
present the methodology and parameters employed to deter- 
mine the anisotropy energy spectrum and the EM cascades. 
In Section 3, we present the results of the calculation. 



2 METHODOLOGY 

2.1 Collective Intensity of Unresolved 
Astrophysical Sources 

The contribution to the EGB due to unresolved blazars can 
be viewed as the superposition of the collective intensity of 
intrinsic blazar spectra and the contribution from the cas- 
cade radiation from the interactions of VHE photons from 
blazars with the EBL: 



Ie(E ) 



Ie(Eo) + Ie(Eo) 



(1) 



where the intensity, Ie(Eo), is given in units of photons per 
unit area per unit time per unit solid angle per unit energy 
emitted at observer frame energy Eq. In calculating the col- 
lective spectrum of unresolved blazars, we follow the formal- 
ism outlined in Venters & Pavlidou (2011). Blazar spectra 
are taken to be smoothly broken power laws: 



Fe{Eo) = Fo 



(I)' 



(2) 



where Fe{Eq) is the differential photon flux in units of pho- 
tons per unit area per unit energy per unit time, E b is 
the break energy, ct\ = a — Aai is the low-energy slope, 
Q2 = ol + Ae*2 is the high-energy slope, a is the gamma-ray 
photon spectral index of the blazar assuming an unbroken 
power-law spectrum, 2 and n quantifies the sharpness of the 
transition from the low-energy power law to the high-energy 
power law (taken to be 1). 

The total flux, F, of photons with energies greater than 
some fiducial energy, Ef (taken to be 100 MeV), is found by 
integrating Fe{Eq) over energy, 



/>oo 

= F 

J E f 



dE . 



(3) 



Then, the contribution of a single unresolved blazar to the 
EGB is 



h = 



F[(E /E b ) a i +(E /E b )^]- 1 
4tt/~ [(Eo/Et)*! + (Eo/Et)^}- 1 dE 



(4) 



where the flux of one source is uniformly distributed over 
the entire sky in anticipation of an isotropically distributed 



cosmological population. Including absorption due to the 
EBL, the contribution becomes 



h = 



FjjEo/EbT + (Eo/E b ) a *}- 1 
**Ib, l(Eo/E b )«i + (Eo/Ef,)^}- 1 dE ' 



-t(E ,z) 



(5) 



where t(E ,z) is the optical depth for a photon with ob- 
server frame energy Eo, emitted at redshift z. In determin- 
ing the total unresolved blazar contribution including EBL 
absorption, one would generally rewrite Equation 5 in terms 
of the blazar gamma-ray luminosity and redshift, multiply 
by the blazar GLF and intrinsic spectral index distribution 
(corrected for the spectral bias endemic to flux-limited sur- 
veys 3 as per Venters et al. (2009)), and finally integrate over 
gamma-ray luminosity, redshift, and spectral index, thereby 
accounting for evolution in both the sources and the EBL 
(through the redshift-dependent optical depth; see e.g., Ven- 
ters et al. 2009). However, in performing a calculation that 
includes cascades, the population of blazars must be divided 
into its subpopulations of FSRQs and BL Lacs since these 
subpopulations have distinct spectral properties. In order 
to perform a full calculation that accounts for the evolution 
of sources (as in Venters (2010)), we would require reliable 
GLFs for each subpopulation. In the case of FSRQs, a reli- 
able GLF does exist (Ajello et al. 2012), but one has yet to 
be determined for BL Lacs since it is difficult to obtain reli- 
able redshifts for them. Given that the purpose of this paper 
is to serve as a proof-of-concept for using relative changes in 
the EGB anisotropy to constrain the IGMF, for the sake of 
simplicity, we do not incorporate full blazar GLFs in our cal- 
culations. Rather, we assume typical luminosities and red- 
shifts for all blazars, deferring a more detailed calculation to 
the future, at which time reliable GLFs for both FSRQs and 
BL Lacs will have been measured. As such, the only integral 
that must be performed is that over the spectral index: 



Ie(Eo) 



-I. 



Io / dap(a) 



[(E /E b )^ + (E /E b )^Y 
S(E f ,a) 



r(.E ,z) 



(6) 



where Iq is a normalization constant, p(a) is the measured 
spectral index distribution (SID), and 



S(E f ,a) 



poo 

= / dE 

Je } 



(7) 



We fix the normalization constant so that the total astro- 
physical (blazars, starforming galaxies, and EM cascades) 
contribution to the EGB fits the spectrum of the EGB as 
measured by the Fermi-L AT (Abdo et al. 2010c). 

In determining the measured SID, we make use of the 
results presented in Venters & Pavlidou (2011) from a likeli- 
hood analysis fitting Fermi-LAT First Catalog (Abdo et al. 
2010e) FSRQs and BL Lacs to Gaussian SIDs accounting for 
errors in measurement of individual blazar spectral indices. 
For these blazar subpopulations, the maximum-likelihood 
Gaussian SIDs can be characterized by means (ao) and 
spreads (ao) with parameters determined to be ao = 2.45 



2 For the sake of clarity, we use a for the photon index rather 
than r, which is commonly used in the literature. 



3 We note that the Fermi survey is not exactly a flux-limited 
survey due to the non-uniformity of the total diffuse background 
throughout the sky. However, for the purposes of this paper, we 
neglect this effect. 
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and cro = 0.16 for FSRQs and a = 2.17 and a = 0.23 for 
BL Lacs. Following Venters & Pavlidou (2011), we model the 
spectral breaks by taking Acn = 0.1, Act2 = 0.9, Et,a = 4 
GeV for FSRQs, and E bfi = 10 TeV for BL Lacs. 

In principle, one would prefer to perform the above pro- 
cedure on subsamples of local universe blazars to ensure that 
only intrinsic spectra are used (i.e., that the spectra have 
not already been modified by EBL absorption and EM cas- 
cades). However, the population of local universe sources is 
sparse, and the results of a likelihood analysis performed 
on such a population would be subject to a considerable 
degree of uncertainty. Even so, in the case of FSRQs, the 
spectra are too soft to produce an appreciable amount cas- 
cade radiation (Venters 2010). For this reason and because 
the photon index of the spectrum of cascades is relatively 
hard (r ~ 2; see Venters 2010 and Section 3 of this paper), 
we do not expect substantial modifications of FSRQ spectra 
due to EM cascades (conversely, we do not expect a blazar 
with a hard measured spectrum to have a much softer intrin- 
sic spectrum). Notably, the spectra of local universe FSRQs 
appear to be similar to that of the FSRQ population at 
large (Abdo et al. 2010e). On the other hand, the spectra of 
BL Lacs are sufficiently hard that one might be concerned 
that the spectra are substantially modified by EM cascades. 
Nevertheless, we note that the population of local universe 
BL Lacs appear, on average, to have harder spectra than 
the BL Lac population at large. As such, rather than at- 
tempting a likelihood analysis on the sparse local universe 
BL Lac population (which also have large measurement un- 
certainties in their spectral indices), we make use of the SID 
for the whole BL Lac population as this represents a more 
conservative choice. 

Since the AGN population peaks ~ 1, we take the typ- 
ical redshift of blazars to be z BIC = 1. In modeling the sup- 
pression of the source brightness due to magnetic deflection 
of cascades (see Appendix), we take the jet opening angle of 
blazars to be ~ 1°. The model of the EBL was taken from 
Franceschini et al. (2008). 

We also include a model of the starforming galaxy (SF) 

I 



where dNr/dE is the spectrum of Inverse Compton scat- 
tered radiation per electron of Lorentz factor, T, P(f; E p , z) 
is the probability that the pair production interaction at a 
given redshift, z, of a primary photon of energy E p will pro- 
duce electron-type particles of energies E e \ — f x E p and 
E e2 = (l-/)x£p,ri = fxE p /mc 2 ,Y 2 = (1-/) x E p /mc 2 , 
and d 2 N~f / dzdE p is the continuous spectrum of photons un- 
dergoing pair production interactions as a function of red- 
shift and primary energy. 

The full numerical integration of Equation 9 is compu- 
tationally cumbersome. As such, in order to calculate the 

I 



contribution to the EGB. In this case, we assume that the 
gamma-ray luminosity of an SF galaxy is mainly due to 
the decay of neutral pions created through galactic cosmic- 
ray interactions with neutral hydrogen in the galaxy. It is 
thought that galactic cosmic rays are accelerated in su- 
pernova remnants and then diffuse throughout the galaxy. 
Hence, the flux of galactic cosmic rays is expected to be 
proportional to the galaxy's supernova rate, which is, in 
turn, proportional to its star-formation rate. Furthermore, 
the neutral hydrogen gas provides the fuel for forming stars. 
As such, the expectation is that an SF galaxy's gamma- 
ray luminosity can be parameterized in terms of its star- 
formation rate (Pavlidou & Fields 2002; Abdo et al. 2010a; 
Fields et al. 2010; Makiya et al. 2011; Stecker & Venters 
2011; Ackermann et al. 2012b). For the purposes of this 
paper, we model the starforming galaxy contribution using 
the infrared luminosity function model of Stecker & Ven- 
ters (2011). In this model, the star-formation rate of an SF 
galaxy is assumed to be proportional to its infrared lumi- 
nosity; thus, the gamma-ray luminosity of an SF galaxy is 
proportional to a power of its infrared luminosity. In Stecker 
& Venters (2011), the relationship between the gamma-ray 
luminosity of an SF galaxy and its infrared luminosity was 
determined by fitting the measured gamma-ray luminosities 
of Fermi-detected SF galaxies to their infrared luminosities. 
The SF galaxy contribution could then be determined by in- 
tegrating over an infrared luminosity function of SF galaxies, 
which they took from Hopkins et al. (2010). 

2.2 Electromagnetic Cascades from Blazars 

The cascade intensity is given by 



fc dtdAdQdE ' v ; 

where dN c /dE is the spectrum of cascade photons due to 
pair production and Inverse Compton scattering (for full 
equations and derivations, see Venters 2010): 



e- T{E °' z) dfdE p dz , (9) 
I 

cascade spectrum, one typically makes use of a Monte Carlo 
code, such as Cascata (Venters 2010), that propagates pho- 
tons from known spectra and redshift distributions and cal- 
culates the cascade spectrum of each photon. However, for 
the purposes of this paper, we employ a more semi-analytical 
model of the cascade development. Since we have assumed a 
typical redshift for all blazars, we have eliminated the need 
for randomly generated source redshifts. Instead, the spec- 
trum of pair-producing photons, d 2 N 1 /dzdE p , can readily 
be determined from the absorbed collective intensity of un- 
resolved blazars binned in energy: 



dN Tl {E (l + zj) dNr a (Eo(l + z)) 
dE dE 



p- z^o^o = /o p-r w (W)a ;,t ( T t)aarl i 1 - dE ° 
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For each energy bin, we determined a characteristic red- 
shift for the start of the cascades based on the criterion 
that the characteristic redshift is that at which the num- 
ber of photons that have been "absorbed" in propagating 
from the source redshift is half the number of the total pho- 
tons "absorbed" in propagating to the present epoch. Using 
this criterion, the characteristic redshift can be calculated 
numerically from 



t*(£o,z*) = In 



e (Eo.Zerc) _|_ ^ 1 



(11) 



Also, we assume that each secondary electron produced in a 
pair-production interaction carries half of the energy of the 
primary photon. Finally, since the number density of CMB 
photons is much higher than those of the other wavelengths 
that comprise the EBL, we assume that the seed photons for 
the Inverse Compton interactions are CMB photons. In so 
doing, we can neglect Klein-Nishina effects in determining 
the spectrum of scattered radiation. The spectrum of scat- 
tered radiation per electron of Lorentz factor, T is given by 
(Venters 2010) 



dN r 
dei 



d 3 N r 
dtdeide 



-dtde = 



d 6 N T , 
dtdeide 



dt 



dE e 



dE e de , 



(12) 



where T' is the electron Lorentz factor at time, t, e is the 
initial soft photon energy, e\ is the scattered photon energy, 
d 3 N r i /dtdeide is the differential scattered photon spectrum 
given by (Blumenthal & Gould 1970) 



d 3 N r , 3<j t c n(e) / , , 



dtdeide 16r' 4 e 2 



2ei In 



4T' 2 < 



(13) 

n(e)de is the number density of soft photons, and \dE e /dt\ 
is the rate of energy loss of the electron given by 

dE e 



4 r ,2 

— - — = orcl 

dt 3 



£ 



en{e)de , 



(14) 



We note that this formalism assumes continuous energy 
losses for the electrons, which is appropriate for the regime 
that we are considering. Energy losses due to synchrotron 
emission are also included, but arc negligible compared with 
Inverse Compton losses for the magnetic field strengths con- 
sidered. Furthermore, this simplified approach has the ad- 
vantage that recycling cascades (in which Inverse Compton 
photons up-scattered to high energies are also allowed to ini- 
tiate EM cascades) is more feasible than in the full Monte 
Carlo approach. As such, for BL Lacs, we also calculated the 
spectrum of second generation cascades. Notably, the effect 
of EM cascades is to redistribute high-energy radiation to 
lower energies; as such, the high-energy component of sec- 
ond generation cascades necessary to produce higher genera- 
tion cascades is greatly reduced. Thus, the contribution from 
generations beyond the second generation is negligible. For 
FSRQs, the intrinsic spectra of the sources are sufficiently 
soft that even the first generation of cascades makes a small 
contribution to their overall collective intensity. 4 The sup- 
pression of the source brightness due to magnetic deflection 
of cascades (see Appendix) is also calculated, and we as- 
sume a uniform magnetic field that evolves solely due to the 



expansion of the Universe [B(z) ~ Bo(l + z) 2 ; Neronov & 
Semikoz 2009], where B Q ~ 5 x 10~ 14 G. We should also 
note that for the purposes of this paper, we assume that 
the correlation length of the magnetic field is much larger 
than electron cooling distances. Alternative scenarios for the 
magnetic field and source spectra will be explored in future 
publications. 

2.3 Anisotropy Energy Spectrum 

The anisotropy energy spectrum, Ce(E), is defined as the 
intensity fluctuation angular power (in units of sr) at a given 
angular scale, £, as a function of energy. In the simple case 
of a two-component background, the total angular power is 
given by (Siegal-Gaskins & Pavlidou 2009) 

CT = f?C ( e 1} + / a a C< 2) + 2/ 1 / 2 Q (lx2) , (15) 

where is the angular power spectrum of component 

(n), /„ = I n (E) / Itot(E) is the energy-dependent fraction of 
the total emission arising from component (n), and (j^ nxm ^ 
is the cross-correlation term for components (n) and (m). 
As Cf ] is a measure of angular fluctuations in units of the 
mean, it is independent of energy for a single population 
of sources with identical spectra. If the two components are 



uncorrelated, c^ nxm ^ 
plified, 



0, and Equation (15) is further sim- 



(16) 



If additional components of negligible angular power con- 
tribute to the total background signal, then these will not 
result in additional terms in Equation (16), but rather they 
will affect C| ot (E) by changing the fractional contributions, 
/, of each of the components that do contribute to the an- 
gular power. 

In the simple scenarios we explore, we always consider 
two components that do contribute to the total angular 
power: BL Lacs and FSRQs. Because in our adopted model 
the two contributions do not add up to the total EGB sig- 
nal measured by Fermi (Abdo et al. 2010c), we assume that 
the rest of the EGB signal originates from a very \ow-Ce 
population, such as starforming galaxies. For simplicity, we 
take C| SF ' = for this component; hence, their contribution 
only enters Equation (16) through the total background in- 
tensity, Itot(E) (i.e., the denominators of fi and f2). As 
discussed in the previous sections, the energy spectrum of 
the primary emission from FSRQs is very soft at very high 
energies, so the cascade signal from FSRQs is negligible; all 
cascade emission that we consider comes from BL Lac VHE 
photons. 

In the limit of low B, there is little angular spread in the 
cascade photons, so these photons appear to originate from 
the same sources as the primary emission. In this case, the 



4 We note that the small amount of cascade radiation is not sim- 
ply the result of the low break energy chosen for this analysis; 
Venters (2010) modeled FSRQ spectra as unbroken power laws 
up to ~ 100 TeV and also found that the FSRQ cascade contri- 
bution was small. 
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only effect of the cascades is to alter the energy spectra of 
blazars, leaving their anisotropy properties unchanged. As 
such, the fractional components of Equation (16) become 
f hl (E) = [llT(E) + I^r(E)] //tot(B), while cf" remains 
unchanged. 

In the limit of strong B (Bo ~ few x 1CF 14 G), we deter- 
mined the scattering angles of the cascades (see Appendix) 
and found that for observed photon energies less than ~ 100 
GeV, the scattering angles were sufficiently large such that 
within the observed angular extent of the cascade emission, 
the number of sources above the Fermi-LAT sensitivity is 
sufficiently large to introduce source confusion 5 . Thus, we 
take the cascade signal as having angular power, cf asc) = 0. 
Since the cascades now have an angular power that is dif- 
ferent from the blazars and null, they will enter Equation 
(16) in the same manner as the starforming galaxies, only 
through the total background intensity, I to t(E). The numer- 
ators of /i and fa in Equation (16) will then only include 
the intrinsic emission from FSRQs and BL Lacs, respectively 
(i.e., f bl (E) = llf(E)/I tot (E)). We emphasize that it is for 
this reason that the total angular anisotropy energy spectra 
of the two magnetic field scenarios will be different. Even if 
the overall EGB cascade contributions in the two scenarios 
were the same, the total angular anisotropy energy spec- 
tra would be distinct, depending on the prominence of the 
cascades in relation to the other contributions to the EGB. 

Since, as stated in Section 2.1, a reliable BL Lac GLF 
has yet to be constructed, a full calculation of the angu- 
lar power of BL Lacs as per Ando et al. (2007) is not 
yet feasible. Instead, we normalize our anisotropy energy 
spectra such that in both scenarios, C| ot ~ 10~ 5 sr around 
1 GeV, consistent with the Fermi measurement of angu- 
lar anisotropies in the EGB at these energies (Ackermann 
et al. 2012a). This results in a blazar angular power 6 of 

qCFSRQ) = C (BLLac) = Q x 1Q - 5 ^ ^ bg tegtable in 

future publications once a reliable BL Lac GLF is found. The 
measurement at 1 GeV was chosen to balance the competing 
needs of adequate angular resolution and sufficient photon 
statistics, and to minimize possible contamination from any 
remaining galactic foreground emission. The uncertainty at 
1 GeV is relatively small, though more data from Fermi will 
likely improve it. 



3 RESULTS 

In Figure la, we plot the spectra of the modeled contribu- 
tions to the EGB. The "total" contribution is the sum of the 
contributions from blazars, SF galaxies, and one of the cas- 
cade models (the total is similar for both models since the 
largest differences in the cascade contributions in the two 
scenarios occur at energies at which the EGB is dominated 
by blazars and starforming galaxies). The blazar component 
hardens above ~ few GeV due to the increased contribution 
of BL Lacs with respect to FSRQs, as BL Lac spectra remain 
hard while FSRQs soften. The cascade spectrum in the case 



5 Since the blazar contribution to the EGB is the collective emis- 
sion of unresolved blazars, this criterion is quite conservative. 

6 Note that for a Poisson distribution of sources, the angular 
power at a given energy is constant with multipole. 



of the zero IGMF is, as expected, quite flat (dN/dE oc E~ 2 ; 
see e.g., Strong et al. 1974; Venters 2010) demonstrating the 
cascade effect of redistributing radiation from high energies 
to lower energies. The sum of the three components fits quite 
well the spectrum of the EGB as measured by Fermi. No- 
tably, the spectra of the cascades for the two cases of the 
IGMF, while similar at the highest energies, are quite dis- 
similar at the lower energies. This is due to the loss of some 
cascade radiation in the strong IGMF case as significant 
amounts of the cascades are deflected away from the line- 
of-sight of the observer. However, it should be noted that 
given our simple model of the blazar contribution to the the 
EGB, we were only able to account for cascade losses due to 
deflection. That is, in our simple picture, we are calculating 
the amount of radiation that an observer looking head-on at 
a given source would be able to detect. It is likely that just 
as cascades from one aligned source are deflected out of the 
observer's linc-of-sight, cascades from other aligned sources 
and even misaligned sources are deflected into the observer's 
line of sight. Nevertheless, we note that even if the strong 
IGMF cascade radiation were comparable to that of the zero 
IGMF case, the angular power in the two cases would remain 
distinct since the angular properties of the cascades are dif- 
ferent (see Section 2.3). 

In Figure lb, we plot the anisotropy energy spectra of 
the EGB for the two cases of the IGMF. Notably, there is 
an appreciable difference between the two cases starting at 
~ few GeV. This is due to the fact that above these en- 
ergies, the contribution from EM cascades relative to the 
other components is much more significant than at lower 
energies. As a result, the impact of the isotropization of the 
cascades by the IGMF is much more significant at these 
energies. Aiding in the distinction is that the cascades aug- 
ment the anisotropy in the case of zero IGMF, because they 
enhance the contribution of the most anisotropic compo- 
nent. For comparison, the Fermi-LAT measurement of the 
anisotropy energy spectrum of the EGB (Ackermann et al. 
2012a) is also plotted. With just two years of data, the Fermi 
measurements seem to already favor the higher IGMF case 
over the null IGMF case. 



4 CONCLUSIONS & DISCUSSION 

We have studied the effect of a strong IGMF on EM cas- 
cades from blazars and the resulting impact on the cascade 
component of the EGB and the anisotropy energy spectrum 
of the EGB. We have shown that while the spectrum of 
cascade radiation in the case of zero IGMF is flat in E 2 , 
the spectrum in the case of a strong IGMF is suppressed at 
low energies, though there could be an additional contribu- 
tion arising from cascades from other sources, such as radio 
galaxies and propagating UHECRs. 

For the two extreme cases of IGMF strength we have 
considered here, we have calculated the anisotropy energy 
spectrum of the EGB. In the case of zero IGMF, the ef- 
fect of cascades is to augment the anisotropy arising from 
blazars since the cascades increase the contribution to the 
EGB arising from blazars and are highly collimated with 
the primary emission. In the case of the strong IGMF, the 
cascades are deflected substantially and the angular proper- 
ties of the parent blazar population are blurred. Thus, in the 
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Figure 1. a: The spectrum of the modeled EGB (thick solid line) together with individual components: intrinsic blazar emission, 
including both FSRQs and BL Lacs (thin solid line); SF galaxies (dashed line); electromagnetic cascades for zero IGMF (dot-dashed line) 
and nonzero IGMF (double dot-dashed). For reference, the spectra of the EGB based on both Fermi (filled circles; Abdo et al. 2010c) 
and EGRET data (open squares; Sreekumar et al. 1998; open triangles; Strong et al. 2004) are also plotted, b: The anisotropy energy 
spectra of the modeled contributions to the EGB. The solid line is the total model assuming zero IGMF. The double dot-dashed line is 
the total model assuming nonzero IGMF. Grey boxes represent the Ferroi-LAT measurement of the anisotropy energy spectrum of the 
EGB (Ackermann et al. 2012a). 



case of a strong IGMF, the effect of cascades is to reduce the 
anisotropy of the EGB as they hinder the anisotropy arising 
from blazars rather than reinforcing it. At energies at which 
the relative contribution to the EGB arising from EM cas- 
cades is significant, the anisotropy energy spectra of the two 
cases are distinguishable, and current Fermi data already 
seem to favor non- negligible IGMF values. With more data 
and further refinement of this technique, Fermi will be able 
to provide tighter constraints and possibly even measure the 
IGMF. 

Regarding the stability of EM cascades during prop- 
agation, Broderick et al. (2012) present the very interest- 
ing notion that plasma interactions between cascades and 
the intergalactic medium would introduce energy losses that 
are more significant than Inverse Compton energy losses, 
thereby quenching the cascades. Unfortunately, this idea was 
premised on the assertion that cascades constitute plasma 
beams, which Broderick et al. suggest as a result of the ob- 
servation that there are a large number of pairs in a given 
volume scale. However, a large number of particles in a small 
scale is not sufficient to constitute a plasma - the particles 
must be able to influence each other through their elec- 
tric and magnetic fields (for this reason, the particles in 
a plasma must necessarily be charged particles). From our 
understanding of relativistic beams of particles, the elec- 
tromagnetic interactions between particles are largely sup- 
pressed (Chao & Tigner 1999). We can see this by consider- 
ing the field around a single charged particle at rest and then 



boosting in the direction of propagation. For large Lorentz 
factors (such as those involved in the EM cascade calcu- 
lations), length contraction compresses the electric and as- 
sociated magnetic fields so that they exist only in a plane 
perpendicular to the direction of propagation of the parti- 
cle. Since the bulk of the particles in the cascades will not 
lie in this plane, they will not be electromagnetically con- 
nected to the particle in question. For the few particles that 
could lie in this plane, the Lorentz forces exerted on them 
largely cancel because they are relativistic and propagate 
co- linearly resulting in near cancelation between the electric 
and magnetic field components. For simple charge distribu- 
tions (e.g., the two-dimensional Gaussian distribution), ana- 
lytic expressions have been derived for the Lorentz forces on 
particles in a beam due to other particles within the beam, 
and such forces have been demonstrated to be suppressed 
by a factor of j~ 2 (Chao & Tigner 1999). As such, even 
for electrons with energies ~ few MeV, the strengths of the 
forces exerted on them are suppressed to the percent level; 
at energies ~ hundreds of MeV and above, these forces are 
negligible. In the case of the nonzero IGMF, the compression 
of the electric and magnetic fields to a plane still applies, and 
because of magnetic deflection, the density of charges falls, 
suppressing the likelihood of finding nearby charges in the 
plane. Therefore, we conclude that the notion of a plasma 
beam does not apply in this case because the particles in 
the cascades cannot exert electromagnetic influence on one 
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another (i.e., cascades do not constitute the necessary col- 
lective phenomena for substantial plasma interactions) . 

Many possible contributions to the EGB have been 
identified in the literature, and a multi-component back- 
ground is plausible in light of current Fermi data. The ex- 
tent to which other contributions will impact the differ- 
ence in anisotropy energy spectra in the two cases consid- 
ered here depends on their angular power and whether they 
also have appreciable cascade contributions. For instance, 
other gamma-ray-emitting radio galaxies for which the jet 
is not aligned with the observer's line-of-sight ( "misaligned" 
blazars) could also contribute to the EGB (Inoue 2011). 
For those radio galaxies that are too misaligned for the ob- 
server to receive substantial intrinsic gamma-ray emission, 
a nonzero IGMF could deflect cascades into the observer's 
field-of-view, giving rise to an additional contribution to the 
EGB, provided that the spectra of these sources are hard 
enough to initiate cascades. These sources will be substan- 
tially fainter than blazars, but are not likely to be as nu- 
merous as starforming galaxies. Hence, they are not likely 
to have substantial angular power and will enter the calcu- 
lation of the total angular power of the EGB in a manner 
similar to that of the starforming galaxies, solely by virtue 
of their relative contribution to the intensity of the EGB. 
We will return to the impact of radio galaxies in a future 
publication. 

In addition to that from radio galaxies, there could also 
be a contribution to the EGB from the annihilation of dark 
matter in the Galaxy. Even if gamma rays from dark matter 
annihilation do not dominate the EGB in any energy band, 
they could have a substantial impact on the total angular 
power of the EGB provided that they make a significant 
contribution to the EGB intensity. The impact of dark mat- 
ter annihilation on both the intensity and angular power 
of the EGB depends on the mass of the dark matter par- 
ticle and the combined cross section for annihilation into 
channels that result in gamma rays. Some models of the 
EGB contribution from dark matter can result in spectra 
that peak at the same energies at which the EM cascades 
contribute most significantly to the EGB. However, the sig- 
nificance of the dark matter contribution to the EGB also 
depends on the amount of substructure present in the dark 
matter halo with the largest signals occurring in scenarios 
in which the halo contains numerous large clumps (Siegal- 
Gaskins 2008). In such scenarios, the angular power of dark 
matter is substantial and will likely enhance the total angu- 
lar power of the EGB at the energies at which the dark mat- 
ter signal peaks (Siegal-Gaskins & Pavlidou 2009). In both 
of the IGMF cases considered in this paper, the impact of a 
substantial contribution from clumpy dark matter at ener- 
gies of tens of GeV would be to increase the angular power 
in the last energy bin in Figure lb. The angular power in 
the two scenarios would still be distinct, though the distinc- 
tion would depend on the significance of the dark matter 
contribution relative to that of blazars and EM cascades. 
Nonetheless, in order to push the two models to yield the 
same angular power in the last energy bin, the contribution 
from dark matter would have to be substantial, and even in 
the case of nonzero IGMF, it could prove difficult to explain 
the flatness in total angular power of the EGB with energy 
as measured by the Fermi-LAT (Ackermann et al. 2012a). 
Even worse is the scenario in which there is no substantial 



cascade contribution, but there is a substantial dark matter 
contribution - without a low-anisotropy contribution to mit- 
igate the total angular power, the angular power in the last 
energy bin would rise substantially. We must also consider 
the possibility that the dark matter halo is smooth rather 
than clumpy. In this case, the angular power of dark matter 
would be minimal, but it could be difficult for dark matter to 
contribute significantly to the EGB. In any case, Fermi has 
already begun to constrain models weakly interacting mas- 
sive particles as dark matter candidates, including tantaliz- 
ing hints of a signal at ~ 130 GeV near the Galactic Center 
(see e.g., Weniger 2012). In the future, the much anticipated 
results from the Large Hadron Collider could provide even 
more clues on the nature of dark matter. Thus, the question 
of the dark matter contribution to the EGB could soon be 
answered. 

Finally, there could also be an appreciable cascade con- 
tribution from the interactions of UHECRs with the EBL. 
The impact of cascades from UHECRs on the total angular 
power of the EGB depends on the still unknown sources of 
UHECRs, but the most likely scenario is that the difference 
in angular power between the two IGMF cases will still be 
present and may, in fact, be augmented. In order for the sep- 
aration in angular power to disappear with the inclusion of 
UHECR cascades, the UHECR sources would have to have 
very little angular power. Such a scenario might have been 
possible in the top-down models of UHECRs, but such mod- 
els are tightly constrained by the upper limits on the UHE 
photon flux determined by the Pierre Auger Observatory, 
and some such models have already been excluded (Pierre 
Auger Collaboration 2009). As such, the current discussion 
of the origins of UHECRs largely plays out in the astro- 
physical arena, and the classes of astrophysical sources that 
are typically thought to be capable of accelerating UHE- 
CRs tend to be rare and powerful sources (as are blazars). 
Hence, the scenario in which UHECR sources exhibit very 
little angular power is unlikely. Furthermore, for electrons 
with energies greater than those considered here, the In- 
verse Compton process becomes more inefficient as Klein- 
Nishina effects become important. Therefore, the cascade 
photons from UHECR interactions that would contribute to 
the EGB are likely produced by electrons with energies simi- 
lar to those found in the VHE gamma-ray cascades produced 
here. Furthermore, UHECRs are themselves charged parti- 
cles and would have already been deflected by the IGMF 
prior to initiating cascades. Thus, the inclusion of UHECR 
cascades would likely also result in a difference in angular 
power between the two IGMF cases, which may even be 
more prominent than demonstrated in this paper. Given the 
possibility of a large contribution to the EGB from UHECR 
cascades (e.g., Ahlers et al. 2010) and the possible implica- 
tions for the sources of UHECRs, we will devote a future 
publication to the study of the impact of UHECR cascades 
on the angular power of the EGB. 

We note that an advantage of our approach is that 
rather than focusing on a few select sources, our method em- 
ploys information from the whole observable sky. As such, 
our approach could prove capable of providing more sensitive 
constraints on the global properties of the IGMF. Further- 
more, our approach is not limited to uniform magnetic fields 
with large coherent lengths; other field configurations can be 
probed with this technique. In practice, in order to use this 
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technique to provide an IGMF measurement, intermediate 
cases of IGMF strength need to be considered as well, and 
the uncertainties entering the calculation of the expected 
anisotropy energy spectra in each case carefully evaluated. 
In the case of intermediate IGMF strengths, a numerical 
simulation of the expected angular power is appropriate as 
the angular anisotropy properties of the cascade radiation 
will be altered with respect to the primary BL Lac emis- 
sion, though not completely suppressed as in the case of the 
strong IGMF. In fact, since the cascade radiation is corre- 
lated with the primary blazar emission, the cross-correlation 
term of Equation (15) would have to be calculated. As such, 
for intermediate IGMF strengths, a simple analytical calcu- 
lation of the cascade effect on the total angular power is not 
possible. Concerning uncertainties, these enter our calcula- 
tion primarily through the model assumptions for each EGB 
component (which can, in turn, self-consistently determine 
their anisotropy properties), and through the EBL model. 
We will expand our investigation of the practical usage of 
the anisotropy energy spectrum of the EGB to constrain the 
properties of the IGMF in future publications. 

In the meantime, this simple, proof-of-concept calcula- 
tion has already demonstrated that the anisotropy energy 
spectrum can be a powerful tool in constraining the IGMF. 
With a few more years of data from Fermi, we can begin 
to probe the IGMF and uncover more clues into the origins 
of ultra-high energy cosmic rays and large-scale structure 
formation. 
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APPENDIX A: SUPPRESSION OF THE 
SOURCE BRIGHTNESS DUE TO MAGNETIC 
DEFLECTION OF EM CASCADES 

In calculating the suppression of the source brightness due 
to magnetic deflection of EM cascades, we must determine 
the angular spread of the cascades, 6 Brc , as viewed by a fun- 
damental observer at the source redshift, z Brc (Figure Ala). 
In order to determine 6 Brc for multiple cascades, we begin 
with considering the problem for one generation of cascades 
(Figure Alb). In this case, 9i, the angular spread of the 
cascade as viewed by a fundamental observer at zi, can be 
expressed in terms of quantities evaluated at the redshift 
of interaction of the "primary" photon, 22, and the angular 
diameter distance from z src = zi to 22, Da {zi, 22): 

0i = n f 2 > , (Al) 
D A (z 1 ,z 2 ) 

where d 2 is the proper lateral distance the electron has trav- 
eled as measured by a fundamental observer at 22 (see e.g., 
Longair 2008), and we neglect Klein-Nishina effects, so the 
electron cooling distance, d 2 , is not cosmological. By geom- 
etry, d 2 is given by 



d 2 = l 2 sin 82 = 2r 2 sin 82 , 



(A2) 



where l 2 is the straight path of d 2 , r 2 is the Larmor radius 
at 22 of an electron with energy, E 2 , in a magnetic field of 
strength, B(z 2 ) = Bo(l + z 2 ) 2 , given by (Neronov & Semikoz 
2009) 



L 
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and 82 — (1/2) 9 2 is the angle through which the electron is 
deflected away from the straight-line path of the "primary" 
photon, which is equal to half of the magnetic deflection 
angle of the electron given by (Neronov &: Semikoz 2009) 



9 2 e = : i- = (3xl0- 6 )(l + 22r 
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(A4) 

The angular diameter distance from z\ to 22 can be found 
from the angular diameter distance from 22 to z\ through 
the Reciprocity Theorem: 

I + 21 



Da (21, 22) 



I + 22 



D A (Z2, 2l) , 



(A5) 



where Da(z 2 , 21) for a flat universe is given by 



Da(z2, 2l) = 



1 + 21 



[Dc(zi) - Do (22) 



(A6) 



and D c (z) is the comoving distance at 2. 

Given 6<i in terms of known quantities, we can find BTC 
for two generations of cascades using an equation similar 
to Equation Al. Returning to Figure Ala, we note that in 
this case, the primary photon is emitted at 2 src and inter- 
acts at 21 emitting electron-type particles that do not prop- 
agate cosmological distances. One electron-type particle is 
deflected away from the straight-line path of the primary 
particle through an angle, Si, and up-scatters CMB pho- 
tons. The up-scattered photons then propagate until they 
interact at 22, creating more electron- type particles. Since 
we neglect generations of cascades higher than the second 
generation, 61 is given by the previously developed formal- 
ism. The angle through which the electron at 21 is deflected 
away from the straight-line path of the primary photon, <5i, 
is also given by Equation A4, but for quantities evaluated 
at zi. 

Given di and 81 , the total angular spread of the cascades 
as viewed by a fundamental observer at 21 is 6\ = 0i + Si . 
The lateral spread of the cascade as measured by a funda- 
mental observer at 22, D 2 , is given by Equation Al: 



D 2 = 9\D A {zi,z 2 ). 



(A7) 



Finally, we can again use Equation Al to find the angular 
spread of the cascades as viewed by a fundamental observer 
at the source redshift: 



D 2 



D A {z BIC , 22) 



(A8) 



The flux suppression of the cascade due to magnetic deflec- 
tion can be found from the flux of the cascade in the absence 
of magnetic fields by relating the solid angle of the deflected 
radiation (determined from 6 BIC ), Af2 d , to that of the un- 
affected emission (in this case, the solid angle of the jet), 
Afi: 



F A n- 



(A9) 



We note that this procedure is only valid for non- 
isotropically emitting sources. Conceptually, there should be 
no flux suppression for an isotropic source because the sym- 
metry of the problem demands that whatever emission that 
is deflected out of the observer's cone of sight is replaced by 
emission being deflected in. 

This paper has been typeset from a T^X./ F/TfjX file prepared 
by the author. 
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Probing IGMF with EGB Anisotropy 




(a) (b) 

Figure 1. a: Geometry for the angular deflection of two generations of EM cascades, b: Geometry for the angular deflection of 
generation of EM cascades. 
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